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ABSTRACT 

We model the asymmetry of the KOI-13.01 transit lightcurve assuming a gravity-darkened rapidly-rotating host 
star in order to constrain the system's spin-orbit alignment and transit parameters. We find that our model can 
reproduce the Kepler lightcurve for KOI-13.01 with a sky-projected alignment of A = 23° ±4° and with the star's 
north pole tilted away from the observer by 48° ±4° (assuming = 2.05 Mq). With both these determinations, 
we calculate that the net misalignment between this planet's orbit normal and its star's rotational pole is 56° ±4°. 
Degeneracies in our geometric interpretation also allow a retrograde spin-orbit angle of 124° ±4°. This is the first 
spin-orbit measurement to come from gravity darkening, and is one of only a few measurements of the full (not just 
the sky-projected) spin-orbit misalignment of an extrasolar planet. We also measure accurate transit parameters 
incorporating stellar oblateness and gravity darkening: i?* = 1.756 ± 0.014 Rq, Rp = 1.445 ± 0.016 Rjup, and 
i = 85.9° ± 0.4°. The new lower planetary radius falls within the planetary mass regime for plausible interior 
models for the transiting body. A simple initial calculation shows that KOI-13.01's circular orbit is apparently 
inconsistent with the Kozai mechanism having driven its spin-orbit misalignment; planet-planet scattering and 
stellar spin migration remain viable mechanisms. Future Kepler data will improve the precision of the KOI- 
13.01 transit lightcurve, allowing more precise determination of transit parameters and the opportunity to use 
the Photometric Rossiter-McLaughlin effect to resolve the prograde/retrograde orbit determination degeneracy. 



Subject headings: techniques:photometric — eclipses — Stars:individual:K0I-13 



1. INTRODUCTION 

The Sun's planets all orbit in planes that vary in orbital 
inclination from one another by just ~ 1° — 7° - a situation 
that supports t he nebular hyp othesis for the solar system's 
origin (see, e.g. Lissaueij[l993l ). The orbital planes all differ 
by between ^ 
and by ~ 0.3° 



3° — 7° from the Sun's equatorial plane, 
— 6.3° from the invariable plane (the plane 
normal to the net solar system angular momentum vector). 
Although angular distances are typically provided for these 



differences, we think that they are better shown in two 
dimensions. Figure [T] shows the location of the orbital 
angular momentum vectors for various solar system objects 
in ecliptic coordinates. 

Given its utility in constraining the origin of the solar 
system, measurement of the mutual inclination between 
extrasolar planets' orbits and/or the inclination of extra- 
solar planets' orbit normal relative to their stars' spin axis 
- 'spin-orbit alignment', ip - can shed light on the origin 
of those extrasolar systems as well. Evidence so far shows 
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Fig. 1. — Relative positions of the angular momentum vec- 
tors for the orbits of solar system planets, along with that for 
the Sun's rotation and the solar system invariable plane (plane of 
the net angular momentum), are plotted in J2000 ecliptic coordi- 
nate s. Planetary orbit l ocatio ns calculated from orbital elements 
from lMurrav fc D crmot't (2000), Sun's rotation pole calculated from 
equation given in Scidclmann ( 1992) . and the invariable plane coor- 
dinates are from lCoxl ||2000| ). 



that while the solar system is not unusual, a wide variety 
of planetary system types exists. 

Some systems have planets that are even more nearly 
coaligned with each other than those in the solar system. 
In particular, the Kepler- 11 six-planet syst em has typical 
mean mutual inclinations of just 1° — 2° (jLissauer et al 
201ll ). As evidenced by Figure [1] 1° — 2° mutual inclina- 



tion is significantly smaller than the mutual inclination of 
the solar system's terrestrial planets. Interestingly, how- 
ever, the solar system giant planets have mutual inclina- 
tions of 0.7° — 2.0°, comparable to those in the Kepler- 11 
system. We think that this similarity in mutual inclina- 
tions between the solar system giants and the Kepler- 11 
planets may indicate similar formation and evolution, but 
distinctly different from that of the solar system terrestrial 
planets. Although coplanarity of a system's planets and a 
planet's spin-orbit alignment are distinct measurements, 
they share the property that the proposed mechanisms for 
misaligning one will likely misalign the other as well. 

Radial velocity measurements of planet-hosting stars 
during transits can indicate spin-orbit alignment via the 
Rossitcr-McLaughlin effect. Rossiter-McLaughlin obser- 
vations have revealed that a majority of close-in extra- 
solar planets appear to be spin-o rbit aligned (20 of 31 
measurements of A tabulated by IPerrvman 2011, page 
129, are within 2a of zero). For example, IWinn et al 



()2006l) showed that HD189 733 has a spin-orbit misalign- 
ment of 1.4° ± 1.1°, which iTriaud et al.l (|2009[ ) later re- 
fined to just 0.85° ± 0.30°. HAT-P-l's measured inclina- 
tion of 3.7° ± 2.1° is consistent with spin-orbit alignment 



(| Johnson et al.l 120081 ) . We note that, based on Figure [TJ 
each of these planets is more closely spin-orbit aligned than 
is Jupiter, which has an inclination of 6.1° with respect to 
the Sun's equatorial plane. Many more planets have been 



determined to be spin-orbit aligned to l ower precision (e.g. 



Simpson et al."20ld; iMiller et al.ll2010l : iHellier et aLll201ll: 



Simpson ct al. 201JJ, to name a few recent measurements). 

On the other hand, Rossiter-McLaughlin observations 
have also revealed a significant population of spin-orbit 
misaligned planets. For instance, HAT-P-30b is misaligned 
by 74° ±9 ° (.Johnson et a l. 2011) an d XO-3 is misalig ned by 
70°±15° iHebrard et al.l[2009l ) (see lPerrvmal] (|201lh Table 
6.1 for a summary of spin-orbit alignm ents). Man y planets 
are even retrograde , like HAT-P-l lb (I Winn et al. 20101^ 



Hirano et al.l I2011D. WASP-17b (Anderson et al 



and HAT-P-7b (Narita et al.l 120091: iWinn et al 



I2OI 



2009bD! 



Winn et al.l (|2010al ) noted that these highly spin-orbit- 
misaligned planets occur preferentially around hot stars, 
implying that planet formation and/or evolution differs 
around more massive stars. 

The Rossiter-McLaughlin effect has also been used 
to d etermine the spin-orbit alignment of stellar bina- 
ries (lAlbrecht et al.l l2007l 120091 l20 111 ). In fortuitous 



cases, a star shows spot activity in the area under which 
the planet transits. Those spots' effects on the transit 
lightcu rve can then be used to constrain spin-orbit align- 
ment (ISanchis-Oieda et al. 2011 : Sanchis-Oieda &: Winn 



20111 iNutzman et al.ll201l[ ). 



Barnes (|2009t ) showed that the spin-orbit alignment 



for planets orbiting rapidly-rotating stars can be deter- 
mined from transit photometry alone, t aking advantage o f 
the nonuniformity of those stars' disks (|von Zeipellll924l) . 
The nonuniformity introduces characteristic asymmetries 
into planets' transi t lightc u rves. The photometric tech- 
nique described by iBarnes has an advantage over 
Rossiter-McLaughlin measurements in that it can mea- 
sure both the longitude of the planet's ascending node 
and the stellar obliquity to the plane of the sky. As 
a purely radial measurement, Rossiter-McLaughlin deter- 
mines only A, the projected spin-orbit alignment. The 
Rossiter-McLaughlin-measured spin-orbit alignment dis- 
tribution must be deproj ected in order t o ascertain the 
true distribution (Fabrvckv fc Winn 2009t ). 



Recently, ISzabo et al. ( 2011 ) found an asymmetry in 
the Kepler transit lightcurve of KOI-13.01 that they at- 
tributed to spin-orbit misalignm ent ar ound a fast-rotating 
star, consistent with the Barnesi ( 2009() predictions. In this 
paper, we use this asymmetry to measure the spin-orbit 
misalignment angle ip of the KOI-13.01 system using the 
iBarnes (2009 ) model. The present work, focused on the 
measurement of (/?, leaves to others analysis of the com- 
plete lightcurve and secondary eclipse. In Section |5| we de- 
scribe data reduction of the Kepler public Q2 lightcurve for 
this planet, and we describe the gravity-darkened rapidly- 
rotating star transit lightcurve model and its application 
in Section |3l We then show the results of our fit in Sec- 
tion m and discuss the implications of the measurement in 
Section |5l We give our concluding remarks in Section [S] 
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OBSERVATION 



(which we call 'full polar'). 



The Kepler mission (jBorucki et al.ll2010a[) science team 
first identified a planet candidate around the star with 
designation 994f66 2 in the Kepler Input Catalog (KIC, 
Brown et al. l201l[ ) i n the first two qua rters of science 



data (Ql and Q2) (|Borucki et all l201l[ ). A discovery 



paper confirming the resulting Kepler Object of Inter- 
est ( KOI) KOI-13.01 a s a planet is presently in prepara- 
tion (jRowe et al. 2011 ). The KIC records the star's lo- 
cation as 19.1314111 -1-46.8684°, and provides stellar pa- 
rameters of M, = 1.83Mq, Teff = 8848 K, and a magni- 
tude in the Kepler bandpass of Kp ~ 9.958 ( Borucki et al.l 
20111). Lat e r add itional ground-based spectroscopy by 



Szabo et al 



(|201lh determined = 2.05Mo, T^s = 
8511 ± 400 K, and a Amagnitude betwe en KOI-13 and 
its companion of 0.29. ISzabo et al. (2011) also measured 



wsini = 65 ± 10 km/s for KOI-13. 

We use the Kepler photometry, in which apparent tran- 
sit depth for KOI-13.01 is different in the QO, Ql, and Q2 
Kepler time series. This interquarter variability (as well 
as intraquarter variability within Ql) derives from use of 
a too-small aperture in Ql that did not gather all of the 
flux from this saturat ed bright star. T herefore for our in- 
vestigation we, like Szabo et al.l ( 2011 ). use just the Kepler 
Q2 public domain photometry, downloaded from the Mul- 
timission Archive at STScI (MAST). We use the Q2 short 
cadence data in our analysis. 

In reducing the data flrst we eliminate data with im- 
perfect trend subtractions that were acquired near data 
gaps, leaving 111858 1-minute measurements in the short 
cadence time series. Then we clip out the points whose flux 
varies by more than 5a away from the 6 surrounding points 
- this eliminates 176 bad measurements. With no evidence 
of transit timing variations or other transit- to-transit vari- 
ability , we then fold the data using the iBorucki et al 
( 201 I n period of 1.763589 days. We crop the time series 



around the transit itself so as not to consider the inter- 
esting effects in the h ghtcurve other than the transit {i.e., 
Shporer et al.ll2011bl ). In order to speed computation time 
in fitting, we bin the data into 60-second bins, typically 
coadding 45 points per bin and leaving 254 points in the 
timeseries with characteristic precision of 28 ppm. 

Finally, we subtract off 45% of the median flux in order 
to compensate for light from KOI-1 3's binary compan ion 
polluting the photometry (following ISzabo et al. 2011 ). 



3. MODEL 

We gener ate model lightcurves usin g the algorithm first 
described in ^Bar nes &: Fortnev ( 2003[ ), suitabl y modified 
for gravity-darkened rotating stars by iBarne i (|20091 It 
has four different methods to calculate transit fluxes. In 
order of increasing accuracy and required computation 
time, they are: analytical (after iMandel &: Agol 20021 
with gravity-darkening modifications), image-pixel- based 
cartesian, radial 1-D numerical integration (which we call 
'polar'), and radial/azimuthal 2-D numerical integration 



For this work we also added the following features to 
the code: 



• Fit for tp, X - We added the ability to fit for the 
stellar obliquity, tp, and for the longitude of the as- 
cending node of the planet's orbit, w. In this pa- 
per we express lu in terms of A, the sky-projected 
spin-orbit angle, in order to provide comparability 
with existing measurements of oth er transits from 
the Rossiter-McLaughlin effect {e.g. Perrvmanll201ll 
page 128). Hence we define A to be the longitude 
of the ascending node as measured counterclockwise 
from the ascend ing node of the star 's rotational equa- 
tor. Given the ISzabo et al. (2011) measurement of 
stellar v sin i (which is v cos ip using our parameters), 
we therefore couple the stellar rotation rate $7* to the 
stellar obliquity such that: 



{v sini] 



measured 



i?* cos 1p 



(1) 



Time integration - When fitting data, for each 
model photometric data point we numerically inte- 
grate the flux for ±t/2 around the central time for 
the point, where t is that point's integration time. 
This has little effect on the short cadence data, but 
greatly improves the accuracy of flts using the 30 
minute effecti ve integration time for Kepler^s long 
cadence data ( Kippind 2010l ). 



Bandpass - We modified our algorithm to inte- 
grate flux numerically across the published Kepler 
response functior0. Wavelength effects are incor- 
porated entirely into limb darkening parameters 
for non-rotating stars. But fast-rotators' temper- 
ature heterogeneity introduces variations in transit 
lightcurve shape as a function of wavelength ( Barnes! 
20091 ). In practice, for KOI-13, on which we fix 



the polar temperature to be 8848K (the value from 
the KIC), the effects of the bandpass are nearly 
indistinguishable from a monochromatic model at 
0.5934 ^m. 



Fit for M, - Transit lightcurves for non-rotating 
stars can constrain only the stellar density, = 
^j^3 . They cannot determine the radius i?* and 
mass individually because, as can be derived from 
Kepler's law, ^ = ^rpl ^ , where G is the universal 
gravitational constant, Iduj. is the transit duration, 
and n is the planet's orbital mean motion (which is 
set by the orbital period). Since G, Tdur, and n are all 
constant, there are no independent constraints on M* 
or i?* individually. For a fast-rotating star, however, 
the intensity of the gravity darkening effect depends 
on the local effective gravity. For instance, the ratio 
of the equatorial stellar effective temperature Toq to 



^http: / /keplergo. arc. nasa.gov/CalibrationResponse.shtml 
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the polar temperature Tpoio is 



rp4 

cq 
pole 



Rl 



(2) 



The 



term provides no new constraints, of course. 

But Q depends only on and not on M^, (-^ is 
also a function of 17 - see ISeager fc Huil (|2002l ) for 
how they relate in a point-core model relevant for 
stellar moments of inertia) . So gravity darkening can 
be used to break the stellar density degeneracy. 

In practice, however, this is a weak constraint. A 
stellar mass of = 0.001 Mq can be immedi- 
ately (if unhelpfully) ruled out for KOI- 13 because 
even with zero obliquity {ip = 0°) it would need to 
be rotating faster than its breakup speed to achieve 
vsini = 65 km/s. Conversely, a Af, — 100 Mq 
star would need to rotate with a sufficiently long pe- 
riod Prot* = 77 that its gravity darkening would be 
negligible. We can derive lightcurve constraints on 
then because low-mass stars have exaggerated 
gravity darkening, and high-mass stars have muted 
gravity darkening. Either can then fail to fit the 
real lightcurve if the gravity darkening differs suffi- 
ciently from that of the real star. Highly asymmet- 
ric lightcurves better constrain on the high end, 
while more nearly symmetric lightcurves like that for 
KOI-13.01 constrain the low-mass end more effec- 
tively. 

Lacking empirical constraints, we fix the gravity dark- 
ening parameter /3 when we fit the lightcurves. We use two 
differ ent assupti o ns for the value of /3: the theoretical value 
from lvon Zeioell ( 19241 ) of /3 = 0.25, and th e experimental 



value of /3 = 0.19 from lMonnier et al.l (|2007n . As expected, 
higher values of /3 lead to more intense gravity darkening 
holding all other parameters equal. The ultimate effect on 
the transit lightcurve of changing /3, therefore, is roughly 
similar to that of varying M^, in that both primarily af- 
fect the intensity of gravity darkening. We found, how- 
ever, that the best-fit values varied less than la when us- 
ing /3 = 0.19 rather than /? = 0.25. We therefore use the 
theoretical value oi (3 = 0.25 in our reported fits. More pre- 
cise photometry and a better determination of the star's 
mass could therefore help to constrain (3 empirically in the 
future. 

We fit the model described above to the data usi ng the 



Leave nberg-Marq uardt algorithm described in Press et al 



(j2007[ ) to efficiently zero in on the best parameters. Lack- 
ing appropriate analytical partial derivatives for most pa- 
rame ters in the gravity-darkened case (though see also lPall 
20081 for analytical derivatives in the spherical star case), 
we calculate numerical partial derivatives with respect to 
each parameter at each point in each step. 

To calculate the 1-standard deviation (Icr) errors on 
each measurement, we follow a modifie d version of the set 
of steps specified bv lPress et al.l ()2007L page 815). Instead 
of using the covariance matrix, we explicitly fit the data 



using a set of fixed values for each parameter. We then 
identified the minimum in a function of each pa- 

rameter. We assign the error for that parameter based 
on how far from the minimum value the best-fit 
creases by the appropriate amount given a Icr confidence 
level and our number of parameters of interest (8). This 
Ax^ = 9.304 in our case, since we held M^, fixed when 
fitting other parameters (see Section HJ. 

4. RESULT 

We fit the Kepler Q2 short-cadence KOI-13.01 transit 
lightcurve for 9 different parameters: 

1. the stellar radius, i?*; 

2. the planetary radius, Rp] 

3. the orbital inclination relative to the plane of the sky, 
i, also expressed as the impact parameter, b; 

4. a single limb darkening parameter, ci, equal to the 
sum of the two quadratic lim b darkening parame ters 
such that ci = ui + U2, after iBrown et al.l (|2001l ): 



5. the time of inferior conjunction, Tq, which can be 
shifte d relative to the mid-transit time for oblate 
stars (|Barnesll2009l ): 



6. the out-of-transit star flux, Fq; 

7. the stellar mass, M^,; 

8. the longitude of the planet's ascending node. A, rela- 
tive to that of the ascending node of the stellar equa- 
tor; and 

9. the stellar obliquity, ip, measured as the tilt of the 
stellar north rotational pole toward Earth relative to 
the plane of the sky. 

We show the best-fit values along with their Icr uncertain- 
ties in Table [T] The fit itself along with the data, are 
plotted in Figure [21 We also show in Figure [H the best 



fit ligh tcurve without gravity darkening, after ISzabo et al 
(|201l[ ). in order to draw attention to the necessity of 



the present analysis. The no-gravity-darkening fit fails to 
match the KOI-13.01 transit ingress and egress. It also 
does not adequately model the transit bottom, which for 
KOI-13.01 is darker in the first half of the transit than it is 
in the second half. The gravity-darkened model fits all of 
these characteristics - the ingress, egress, and asymmetric 
transit bottom. 



The best-fit value for stellar mass AI^, is highly uncertain 
due to the weak dependence of gravity darkening on My 
as discussed in Section [3l Both the iBorucki et al.l (j201l[ ) 
and Szabo et al. ( 201 if ) spectroscopic measurements of , 
= 1.83 Mq and = 2.05 Mq respectively, lie within 
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Best-Fit Values 


Pcira-iTieter 


Mh. = 1.83 Mq 


Mh. = 2.05 M© 


•^rcducca 


1.409 


1.419 




1.694 ±0.013 Rq 


1.756 ±0.014 Rq 


Rp 


1.393 ±0.015 Rjup 


1.445 ±0.016 Rjup 




0.084508 


0.084513 


i 


85.9° ± 0.4° 


85.9° ± 0.4° 


b 


0.31962 


0.31598 




0.49 ± 0.03 


0.48 ±0.03 


Tn 


4628032 ± 3 s 


4628033 ± 3 s 


A 


24° ± 4° 


23° ± 4° 




-45° ± 4° 


-48° ± 4° 




54° ± 4° 


56° ± 4° 


p 

rot* 


22.5 hr 


22.0 hr 


/* 


0.018 


0.021 



Table 1: Best-fit transit parameters for the KOI-13.01 system. 
The time of inferio r conjunction, Tp, is measured in seconds after 
B JD 2454900, after iBorucki et al There is a 4-way degen- 

eracy for the transit geometry, i.e. the projected spin-orbit angle A 
and the stellar obliquity ip - see Figure |3] The stellar rotation period 
is denoted Prot*. The stellar dynamical oblateness is denoted /«. 



2a of our measured value of Af, = 0.9 ± 0.6 Mq. Our 
measurement clearly prefers lower values for A/*. However, 
the best-fit value of Af, =0.9 M© is clearly inconsistent 
with the star's spectral type. Given the uncertainty in 
our measurement, however, for the remainder of the fitted 
parameters we s how values for an a ssumed mass of both 
= 1.83 M.7^ (iBorucki et al.l l20n"h and A/, = 2.05 M© 
(ISzabo et al.ll201]E 



The reduced (Xreduced) values for each of those fits 
are x?educed = 1-409 (M,. = 1.83 M©) and x?educcd = 1-419 
{M^ = 2.05 M©). We think that the X?educed above 1.0 
results from inter-transit variability in the Kepler pho- 
tometer. Although each individual data point from the 
Kepler Q2 dataset has high precision, red noise on longer 
timescales leads to systematic errors between transits. 
Hence folding the transits introduced a modest amount 
of additional noise over and above the precision of the in- 
dividual points. 

We adjusted the formal la error bars to account for the 
final Xreduced ol "^ach fit. The errors that we indicate in 
Table [T] represent the formal errors of the fit. They do 
not include, for instance, likely systematic errors from our 
assumptions of stellar mass, Tpoic, the binary flux fraction, 
and/or fixed ?;sini. 

Our measured value for the stellar radius that we de- 
rive from the lightcurve fit, i?* = 1.756 ± 0.014 R©, is 
su bstantially lower th an the spectroscopic value reported 
bv lSzabo et al. ( 2011 ). Given a similar planet-star radius 
ratio, then, we also measure a smaller radius for the tran- 
siting body, Rp = 1.445 ± 0.016 Rjup. This smaller ra- 
dius for the transiter places it within the range of pos- 
sible planetary co mpanions for relati vely young, highly- 
irradiated p lanets (iFortnev et al.ll2007l) . According to Fig- 
ure 5 from iFortnev et all (|2007l ). both Afp = 1 Mj and 
A^p — 3 Mj planets with 25 M© cores could have radii of 



0.0002 
.0.0001 



Time from Inferior Conjunction (liours) 
-2-10 1 2 




1^^0.0001 

^l|o.oooo 

trtto.OOOl 

-°-0002_2 -10 12 
Time from Inferior Conjunction (hours) 

Fig. 2. — This plot shows the transit lightcurve for KOI-13.01 in 
filled black circles at in the middle, with a best-fit line using gravity 
darkening {Mt, = 2.05 Mq) in blue and one without gravity dark- 
ening in red. The size of the data circles corresponds to their error 
(as reported in the MAST, adjusted for binning) in the j/-dimension 
(flux). We plot the residuals for the gravity-darkened flt as blue dots 
in the bottom graph. The residuals for the non-gravity-darkened fit 
are at top, as red dots. 



i?D = 1.44 R . TuD wi th an age of 700 Myr as determined by 



Szabo et al.l (|201l[ ). This raises the possibility that KOI- 
13.01's mass lies in the pla netary regime , not t hat of brown 
dwarfs as the higher Rp of Szabo et al.l ( 201ll) suggested. 



By fitting for both i?, and our limb darkening parameter 
ci, we arrive at a more consistent measure of the planet's 
orbit inclination relative to the plane of the sky, i. Ex- 
pressed as an impact parameter 6, our values near b — 0.32 
are significantly lower (i.e. the transit is more nearly cen- 
tra l, and less cl o se to grazing) than the value determined 
by ISzabo et a l.' ('2011^. The discrepancy may lie in the 



assumption by Szabo e t al.l (|201lh of a spectroscopically- 
determined stellar radius that was rather higher than our 
best- fit value. Because the stellar oblateness, affects 
ingress and egress ti mes (similar to the way th at planetary 



oblateness does, see iBarnes fc Fortnev our mea- 



sured impact parameter is more accurate than one that 
would result from a fit without explicitly incorporating 
stellar oblateness. 
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Stellar oblateness /* also affects the measured time 
of inferior conjunction, Tq. For rapidly-rotating, oblate 
stars, the mid-transit time and the time of inferior con- 
junction can differ (Barne'il2009l) . depending on the tran- 
sit geometry. The measured times for the fit assuming 
Mh, = 1.83 Mq (4628032 ± 3 s after BJD 2454900) and 
that while assuming = 2.05 Mg (4628033 ± 3 s after 
BJD 2454900) are slightly different because the intensity of 
the effect depends on the stellar oblateness, which is differ- 
ent in the two fits (see below) . Both values are somewhat 



later than the mid-transit time reported in iBorucki et al 
(|201ll ). 4628014±10 s. Since the values are within 2a, how- 



ever, the discrepancy could be due to the stellar oblateness, 
but may also be statistical variation. 

The transit geometry determinations represent the most 
important measurements that we describe in this paper. 
Our value for the sky-projected spin-orbit angle A can be 
directly compared to measurements of A made in other 
transiting systems from the Rossitcr-McLaughlin effect. 
Unlike Rossiter-McLaughlin measurements, however, we 
are also able to determine the stellar obliquity, V'j which 
we define as the amount by which the star's north pole is 
tipped out of the plane of the sky toward the observer. 

We cannot determine A and -0 uniquely, though. We 
determine 4 allowed transit geometries, as diagrammed in 
Figure[31 For the Af* = 2.05 Mq case, those are (clockwise 
from lower-left): (1) A = 24°, tp = —48°, the prograde 
case with north pole tipped away from the observer; (2) 
A = —24°, tfj = 48°, the prograde case with the north 
pole tipped toward the observer; (3) A = —156°, -0 = 48°, 
the retrograde case with the north pole tipped toward the 
observer; and (4) A = 156°, tp = —48°, the retrograde case 
with the north pole tipped away from the observer. 

The combination of A and tp, along with the planet's 
orbit inclination i, allow us to calculate the complete rela- 
tive angle between the stellar rotation pole and the plan- 
etary orbit norma l - the spin-orbit alignment, if - using 
Winn et al.l (j2007f ) Equation 7. With our angle definitions, 
this equation becomes 



cos (f — sin Tp cos i + cos tp sin i cos A 



(3) 



We measure if = 54° ± 4° for M, 
(f = 56° ± 4° for = 2.05 Mq. 



= 1.83 Mq, and 
Hence the orbiting 
body, KOI-13.01, is far from spin-orbit alignment. While 
the projected alignment A (as would be measured using the 
Rossiter-McLaughlin effect) alone reveals a significant de- 
viation from alignment, the full extent of the misalignment 
manifests from the large stellar obliquity tp which would 
not be measured by the Rossiter-McLaughlin effect (i.e., 
Fabrvckv fc Winnll2009l ). 



In actuality, given the orbital geometry degeneracies 
described by Figure [3l (p also has another solution - 
one where the planet orbits retrograde. In that case, 
ip = 126° ± 4° for Ah = 1.83 Mq, and = 124° ± 4° for 
Af* = 2.05 Mq. We offer a possibility for resolving the pro- 
grade/retrograde degeneracy photometrically in Section [5] 



-2 

O 



o 
o 



Prograde 




Retrograde 



J.=-24° 
1^=48' 




A=204' 
>|;=48 





<|/=-4!}' 



Fig. 3. — Allowed geometries for the KOI-13.01 transit event. Best- 
fit values for an assumed M, = 2.05 Mq are shown. We enhance the 
contrast of the stellar image in order to make the heterogeneity of 
the disk more evidently visible. Note that due to the competing 
effects of limb darkening and gravity darkening, the stellar poles are 
not the brightest part of the stellar disk. Gravity darkening does, 
however, displace the brightest part of the disk poleward from the 
center of the disk. "NP" and "SP" stand for North Pole and South 
Pole, respectively. 



5. DISCUSSION 
5.1. Implications 

Our determination of a high spin-orbit misalignment 
for KOI-13.01, ^ 55° rather than the solar system's more 
typical < 6°, drives us to pose the question of how the 
spin-orbit misalignment came to be. Assuming that KOI- 
13.01 formed like a solar system planet, in a planar disk 
orbiting the protostellar KOI-13, there are two ways to 
have developed a spin-orbit misalignment: either the stel- 
lar spin axis changed, or the planet's orbit did (or both). 
In addition to explaining the the spin-orbit misalignment 
of KOI-13.01, any explanation for its formation and evo- 
lution must also explain KOI-13. Ol's apparently circular 
orbit, as establish ed by the presenc e of a properly-timed 
secondary eclipse (jSzabo et al.ll201ll ). 

iBate et aL ( 2010() suggest a mechanism by which to al- 
ter the star's rotation axis, moving it away from the normal 
to the protoplanetary disk. This mechanism could drive 
spin-orbit misalignment, and allows for KOI-13.01 to end 
up on a circular orbit. 



Two prominent ways to alter the planet's orbit are by 
planet-planet scattering and the Kozai mechanism. If a 
planet-planet scattering event occurred early enough in the 
planet formation process, then sufficient debris might have 
remained to circularize the orbit before the disk dissipated. 

A Kozai-driven sour ce for spin-orbit misalignment, 
like that for HD80606b (|Wu fc Murravll2003t I Winn et al 



2009ar ). is appealing because of the presence of KOI-13's 
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binary companion. However this mechanism requires that 
the planet's orbital precession due to the binary compan- 
ion be mor e rapid than that indu ced by the planet's parent 
star alone. IWu fc Murray! (|2003l ) suggest that HD80606b's 
Kozai variations thus ended when its orbital semimajor 
axis dropped, reducing Kozai's effectiveness while increas- 
ing the effects from its parent star. 

For this mechanism to be viable for the evolution of 
KOI-13.01, its orbit would have had to have been cir- 
cularized from a presumed prior highly-eccentric Kozai- 
derived orbit. Given KOI-13's young 700 Myr age, though, 
tidal circularization would seem difficult. The conventional 
tidal circularization timescale Tc = g^ygj ([Trilling 2000j) 
for KOI-13.01 in its present state is Tc =^ 8 Myr. How- 
ever in a presumed former state like that of HD80606b 
presently, with a semimajor axis of 0.5 AU, the circulariza- 
tion time scale is Tc = 200 trillio n years using the value for 
Mp from lShporer et all (2011b). Hence we think that the 
Kozai mechanism is likely not responsible for the spin-orbit 
misalignment of KOI-13.01, since the end state of Kozai 
migration would be a highly eccentric orbit that could not 
have been circularized in the star's lifetime. A more com- 
plete numerical study of the orbital circularization of KOI- 
13.01, tracking the coupled evolut ion of bot h eccentricity 
and semimajor axis following Jackson et al.l (|2008.') , might 
shed additional light on the issue. 

5.2. Future Work 

The present work used just the public Q2 short-cadence 
Kepler data for KOI-13. Presuming that Kepler continues 
to observe KOI-13 at short cadence, the ultimate coadded 
photometric precision should continue to improve. The re- 
sulting higher signal-to-noise ratio will allow more precise 
determinations of all of the transit parameters in Table 
[TJ As the precise shape of the asymmetry at the bottom 
of the lightcurve becomes clearer, the purely photometric 
determination of the star's mass Af* should approve com- 
mensurately. 

Another effect that will become important as precision 
i mproves is the Pho t ometric Ros siter-McLaughlin effect 
(jShporer et al.ll2011at lGrootil2011l ) - the photometric ana- 
log of the Rossiter-McLaughlin effect. This effect arises 
due to the rotation of the star combined with the beaming 
effect. During transit the planet covers up varying parts 
of the star, with different apparent radial velocities. The 
different radial velocities cause photons emitted from dif- 
ferent parts of the stellar surface to be beamed by a differ- 
ent amount, resulting in an anomalous photometric signal 
during transit. This combined effect has been proposed 
to investigate th e spin-orbit alignmen t of Kepler'' s binary 
star population (jShporer et al.ll2011al ). 



We show a calculation of the intensity of this effect for 
KOI-13.01 in Figure |4l While its present amplitude of 
^ 4 ppm renders it undetectable in the present Q2 Ke- 
pler data, future improvements in precision as the mis- 
sion progresses will make the effect detectable and more 
easily discernible from the gravity-darkening effects ad- 



dressed in this paper. Most importantly, the sign of the 
Photometric Rossiter-McLaughlin effect would be flipped 
for retrograde-orbiting planets. Hence the Photometric 
Rossiter-McLaughlin effect holds promise in the future for 
a purely photometric resolution to the prograde/retrograde 
degeneracy that we show in Figure |3l 




ffl 



-1 1 

Time from mid eclipse [liour] 

Fig. 4. — This plot shows the intensity of the Photometric Rossiter- 
McLaughlin effect for KOI-13.01 using our best-fit parameters for 
Mf = 2.05 Mq from Table [T] While our algorithm presently does 
not account for this effect, it could be incorporated in the future as 
additional Kepler data improves the period-folded photometric pre- 
cision. If it could be detected, the Photometric Rossiter-McLaughlin 
effect could be used to break the prograde/retrograde degeneracy for 
the planet's orbit depicted in Figure|3l 



Finally, iBarnes (l2009l) showed that the lightcurve shape 
varies for spin-orbit misaligned planets when viewed at dif- 
fering wavelengths. We show an illustration of this effect 
for KOI-13.01 in Figure [SJ Because the intensity of the 
difference between the lightcurve in the Kepler bandpass 
and that at different wavelengths depends on Tpoio, such 
observations could serve to place constraints on Tpoio. 

6. CONCLUSION 



Earned (|2009l) first discussed the prospect for spin-orbit 
misaligned planets orbiting rapidly-rotating stars having 
asymmetri c transit lightcur v es owing to stellar gravity 
darkening ( von Zeipeilll924 ). ISzabo et al (I2OIII) showed 
that pla net candidate KOI-13 .01, first identified in Kepler 
data by iBqrucki et aD (l201ll) . has an asymmetric transit 
lightcurve. ISzabo et al.l (120111 ) also made new critical mea- 
surements of KOI-13, including the binary flux fraction and 
V sin i from stellar rotation. 

Fitting the Q2 s hort cadence Kepler photometry for 
KOI-13.01 with the iBarned (|2009l) model, we measured a 
projected spin-orbit alignment A = 23° ± 4° and a stellar 
obliquity (defined as the angle by which the stellar north 
pole is tilted out of the plane of the sky toward the ob- 
server) of -0 = —48° ± 4°. While thes e measurements as - 
sume = 2.05 M© as determined bv lSzabo et al. ( 2011 ). 
other fits using different assumed stellar masses show that 
the overall spin-orbit misalignment derived for the star, 
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Fig. 5. — Here we show the difference between various theoreti- 
cal lightcurves, conforming to the transit parameters from Table [T] 
acquired at various indicated monochromatic wavelengths relative 
to that integrated across the Kepler bandpass. At shorter wave- 
lengths the KOI-13.01 transit lightcurve asymmetry is increased, 
and at longer wavelengths the asymmetry is more muted. A single- 
wavelength lightcurve at 0.5934 fim most closely replicates that of 
the lightcurve under the entire Kepler bandpass for KOI-13. Obser- 
vation of this transit, if sufficiently precise, could constrain Tpoic and 
the gravity-darkening parameter /3. 



(fi = 56° ±4°, is not substantially affected by plausible 
stellar mass variations. The photometric determination 
also allows for a retrograde orbit with ip — 124° ± 4°. Our 
measurement of ip is the first spin-orbit determination to 
come from gravity darkening. It is one of just a handful of 
observations of any kind to determine the complete spin- 
orbit misalignment of a transiting planetary -class body, 



not iust its projected spin-orbit al i gnment A (IWinn et al 



til 

20071: ISanchis-Oieda fc Winnll201lUNutzman et al.ll201l[ ). 



The spin-orbit misalignment of KOI-13.01 follows the 
trend that planets around more ma ssive stars are mor e 



likely to be spin-orbit misaligned (jWinn et al.l l2010a[) . 
Our KOI-13.01 measurement together with spin-orbit 
misaligned WASP-33b aroun d a M* = 1. 5 star 



Collier Cameron et al. 2010[ ) extends the Winn et al 



2010aj ) to yet-higher stellar masses. 



We also measure the other critical transit parameters 
for the KOI-13.01 system. We determined i?* = 1.756 ± 
0.014 Rq, Rp = 1.445 ± 0.016 Rjup, i = 85.9° ± 0.4°, 
b = 0.31598, and a single limb darkening parameter of ci = 
0.48 ± 0.03. Our photometrically-determine d stellar and 



planet ary radii are much lower than those of iSzabo et al 
(1201 ih who used a spectroscopically-derived star radius 
of 2.44 Rq. Our lower estimated radius places KOI- 
13.01 in the regime where interior models allow for high ly- 
irradiated planetary- mass objects ( Fortnev et al. 20071) . 



darkening on the stellar surface gravity to directly measure 
A/* from the Kepler transit lightcurve. Transit lightcurves 
for spherical stars can only constrain the stellar density 
p*, and do not provide any absolute length scale with 
which to constrain M*. Smaller stars would need to ro- 
tate faster to provide the measured tisini = 65 km/s, and 
thus they drive more intense curvature at the bottom of the 
lightcurves of transiting objects. Conversely, more massive 
stars have correspondingly lower gravity darkening, and 
less lightcurve asymmetry. The effect is weak enough that 
our highly uncertain measureme nt of M^. = 0. 9 ± 0.6 Mq 



We also take advantage of the dependence of gravity 



cannot differentiat e between the iBorucki et al.l (|201ll ) and 
Szabo et"all ( 2011 ) spectroscopic measurements of = 
1.83 Mq and = 2.05 Mq respectively. While it is 
consistent with (or within 2a of) either measurement, and 
inconsistent with the measured spectral type (A), we do 
weakly prefer lower values of M*. Future observations, if 
more precise, can help to more tightly constrain this first- 
ever measurement of Af, from gravity darkening. 

Additional Kepler photometry may also show evi- 
dcnc e for the Photometric Ros siter-McLaughlin effect 
(Sh porer et al.ll2011aHGrootll201lh in the KOI-13.01 tran- 
sit lightcurve. This sign of this effect, if it can be mea- 
sured, would resolve the degeneracy between prograde and 
retrograde planetary orbits using photometry alone. 

We have not determined the cause of KOI-13. Ol's spin- 
orbit misalignment. However we note that the Kozai 
mechanism, implicated for the evolution of IID80606b 
(Wu & Murrav 2003), may be inconsistent with KOI- 
13. Ol's circular orbit according to a simple analysis, de- 
spite the presence of a similar-mass binary companion to 
KOI-13. A more sophisticated study of the effectivity of 
the Kozai mechanism for KOI-13.01 can help to resolve 
the origin of the spin-orbit misalignment. Stellar spin drift 
and planet-planet scattering during formation remain vi- 
able misalignment generation mechanisms. 

The complex and interesting nature of the KOI- 
13.01 star-planet system make it an astrophysical labo- 
ratory for photometric transit effects. The exquisite pho- 
tometric precision of the Kepler photometer allows us to 
probe those effects for the first time. The lessons that 
we learn here will then be applicable to future transiting 
objects discovered by Kepler and others. 
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